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The issue of the thermodynamics of a system of distinguishable particles is discussed in this paper. 
In constructing the statistical mechanics of distinguishable particles from the definition of Boltzmann 
entropy, it is found that the entropy is not extensive. The inextensivity leads to the so-called Gibbs 
paradox in which the mixing entropy of two identical classical gases increases. Lots of literature 
from different points of view were created to resolve the paradox. In this paper, starting from the 
Boltzmann entropy, we present the thermodynamics of the system of distinguishable particles. A 
straightforward way to get the corrected Boltzmann counting is shown. The corrected Boltzmann 
counting factor can be justified in classical statistical mechanics. 
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INTRODUCTION 

The concepts of distinguishable and indistinguishable 
particles is important in Statistical Mechanics as their 
corresponding entropies are different. The entropy in sta- 
tistical mechanics is denned in terms of the logarithm of 
the number of the accessible microstates in the phase 
space. The definition of the entropy is called Boltz- 
mann entropy in which it is adopted in popular textbooks 
[3, 0, 01, El- The microstates numbers for distinguish- 
able and indistinguishable particles are certainly differ- 
ent and then their corresponding Boltzmann entropies 
are different, too. However, it is not trivial to know 
whether the distinguishability property may lead to dif- 
ferent physical results. For example, the Gibbs paradox 
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[H, IMlM [il, (H, |I3 |H [U [U [Hf presents one of the 
cases. 

In its simplest case of the Gibbs paradox, the entropy 
of two identical volume of (distinguishable) ideal gas in- 
creases after mixture in which it violates our intuition. 
Consider two subsystems of equal volumes V and parti- 
cle numbers N are separated by a wall, the Boltzmann 
entropy of one system is 



S = klnV 



N 



(1) 



Hence the total entropy Si of two systems is double, 

Si = 2khiV N (2) 

in which the additivity of the Boltzmann entropy is as- 
sumed. Now the wall is removed, the gases of two systems 
mix. The entropy 5*/ is 



S f = kln(2V) 2N 

Thus the mixing entropy reads 

AS = Sf-Si 
= 2Nkln2 



(3) 



(4) 



which is positive meaning that the mixing process is irre- 
versible according to the second law of thermodynamics. 
Here we also assume the Boltzmann entropy is equiva- 
lent to the usual thermodynamical entropy (that is, the 
entropy identified in thermal properties). 

The essence of the paradox is in fact that the entropy 
of the distinguishable particles is not extensive. The en- 
tropy of the mixture is not equal to the sum of their par- 
tition. To resolve the Gibbs paradox, one introduces the 
indistinguishable particles in which a permutation factor 
1/N\ (N is the total particle number) is included in the 
total microstate number to overcome the overcounting 
[(]{. Hence, the entropy before mixing in Eq.([2]) should 
be modified as 



Si = 2k In 



N 



V 
Nl 



(5) 



and similarly, the entropy after mixing in Eq.([3]) reads 

({2V) 2N 



Sf = kin. 



The mixing entropy is 



V (2 AO! 



(6) 



AS = 2Nk\n2-k\n 



OQnN) 
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(7) 



In the thermodynamical limit N — > oo, the mixing en- 
tropy per particle vanishes, and hence the Gibbs paradox 
is resolved. 

The concept of indistinguishable particles is trivial in 
Quantum Mechanics. The identical particles in Quan- 
tum Mechanics are indistinguishable particles. Strictly 
speaking, one cannot distinguish two identical particles 
after their collision. At equilibrium, the number of mi- 
crostates of the whole system is described in terms of the 
number representation (or the second quantized repre- 
sentation). The resolution of the Gibbs paradox is then 
straightforward in views of the quantum nature of iden- 
tical particles. 
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However, we have also similar case in Classical Me- 
chanics, for example, colloids. The colloidal particles 
(giant molecules) of size up to microns are distinguish- 
able since its behavior should be governed by classical 
mechanics. In such a classical system, the entropy is cer- 
tainly not extensive [2,0] and the Gibbs paradox appears. 
Traditionally, to avoid the Gibbs paradox, colloidal parti- 
cles are treated as indistinguishable particles [33| without 
any explicit reason. There is still a puzzle even though 
this kind of treatment turns out to be correct. 

In this paper, we present the thermodynamics of the 
system of distinguishable particles, starting from the def- 
inition of entropy for the distinguishable particles. The 
presentation can easily show how the "reduced" entropy 
[321 ] instead of the original entropy determines the ther- 
modynamical behaviors. 

Before introducing our treatment for distinguishable 
particles, we classify the system of particles into three 
categories, 

1. Distinguishable particles of the same species 

2. Indistinguishable particles (certainly of the same 
species) 

3. Particles of different species (certainly distinguish- 
able) 

by their phase spaces. Suppose we call the phase space 
of the system of N distinguishable particles of the same 
species (the first category) be Tjv. The phase space of N 
indistinguishable particles (the second category) is then 
^n/Sn with the permutation group Sn- The permu- 
tation group Sn is used to eliminate the overcounting 
microstate numbers from Tn. 

For the third category, each particle belongs to its par- 
ticular species. The corresponding phase space of N par- 
ticles of different species is Ti (g) T± <8> . . . ® Ti (total N 
direct products). 

The thermodynamics for the second and third cate- 
gories were already well discussed in the textbook Q. 



THERMODYNAMICS OF DISTINGUISHABLE 
PARTICLES 

To study the thcrmodynamical variables of the distin- 
guishable particle system, we consider two subsystems of 
particle number Ni, iVjj, volumes Vi, Vz, energies Ei,E 2 , 
respectively, in which their particles, volumes, and ener- 
gies are allowed to be exchanged. 

Before two subsystems in contact, the phase space 
is Tjvi ® rjv 2 - After contact, the phase space be- 
comes TjVi+jVa which is larger than <£> IV 2 by (JVi + 
N 2 )l/(Ni\N 2 \) times. This number is the number of way 
to select N\ particles and N 2 particles from the total 
iVi + N 2 particles. That is, 



in which there are total (Ni + N 2 )\/(N 1 \N 2 \) copies 
(or configurational degeneracy). This degeneracy was 
pointed out by Penrose [||. 

The total entropy as a whole becomes 



Si(JVi, V u Ei) + S 2 {N 2 ,V 2 ,E 2 ) + fcln 



(jVj + JVg)! 
N\N 2 \ 

(9) 



Si and S 2 are the entropies of the subsystems separately. 
The last term arises from the configurational degener- 
acy. The existence of the additional term is due to the 
nonextensive property of the system. The distinguisha- 
bility implies nonextensivity of the entropy [3j] in which 
the total system entropy is not the naive addition of the 
subsystem entropies. Similar cases are also well known 
in other physical system [33| . 

Because of that, it is easily seen that the Gibbs paradox 
usually mentioned in the textbook |2j vanishes when we 
look at the phase space of the system of distinguishable 
particles carefully. 

When two subsystems are in equilibrium, the total en- 
tropy attains its maximum, then we have 



S = S 1 (N 1 ,V 1 ,E 1 ) 



S 2 (N 2 ,V 2 ,E 2 ) + S (N U N 2 ) 

(10) 



with So = fcln((iVi+iV 2 )!/(iVi!7V2!)) fromEq. ©. Under 
the constraint of energy E\ + E 2 — E, volume V\ + V 2 = 
V, and particle number N± + N 2 = N, the equilibrium 
attains when the entropy becomes extremum such that 



dS 
~dE~i 
dS 
dVi 
dS 



dSi dS 2 
'dE 1 
dSi 
dVi 
dSx 



dE 2 
dS 2 
dV 2 ' 
dS 2 







dSo_ 

dNi dNx dN 2 ^ dN x 







(11) 
(12) 

(13) 



We can then obtain, for the system of distinguishable 
particles, the expression of the temperature T as 



1 

T 



dS 1 



dE x 

the pressure P as 

P _ f dSi 

By noticing that 

dS Q d 



Ei,Ni 



dS 2 
dEo 



dS 2 
dV 2 



d 



V 2 ,N 2 



E 2 ,N 2 



(14) 



(15) 



dNi 



dNi 



(fclniV 1 !) + — (k\nN 2 \) (16) 
oN 2 



we can express the chemical potential /i as 



T 



diSi-khxN^.) 



dNi 



d(S 2 -k\nN 2 \) 



8N 2 



E 2 ,V 2 



F Nl +N 2 = {T Nl ® rjv 2 ) © (T Nl ® T N2 ) 



(8) 



(17) 
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Now we introduce the "reduced" entropy [3z 



S 1 



S-k\nN\ 



(18) 



such that the above expressions for thermodynamical 
variables 



fdS 



red 



i 



T \ BE 



V,N 



P 

T 

£ 
T 



dS r 



dV 



E.N 



dS Icd 
ON 



(19) 
(20) 
(21) 



E,V 



are re- written in the usual way of thermodynamics for in- 
distinguishable particles. For the system of distinguish- 
able particles, the "reduced" entropy 5' rod instead of the 
original entropy S defines the thermodynamical variables 
of the system. Our presentation is somehow a straight- 
forward way to justify the "reduced" entropy adopted 
for the classical (distinguishable) systems in condensed 
matter physics [341 ]. 



For the case of indistinguishable particles defined in 
the phase space Tn/Sn, its entropy is simply equivalent 
to the "reduced" entropy mentioned above. Although 
the original entropies for both the distinguishable and in- 
distinguishable particles are different, the entropies gov- 
erning their corresponding thermodynamics are still the 
same. Thermodynamics cannot tell the distinguishability 
of the system. 

Suppose our system is N\ , V% , E\ in contact with the 
reservoir of N 2 , V 2l E 2 in which Ni <C jV 2 , V\ < V 2 , E\ < 
E 2 , the entropy of the reservoir can be expanded into 
Taylor's series around N, V, E such that 

S = S 1 (N 1 ,V 1 ,E 1 ) + S 2 (N 2 ,V 2 ,E 2 ) + S (N 1 ,N 2 ) 
= Sl ed (Nx, V 1 ,E 1 ) + S r 2 cd {N 2 , V 2 , E 2 ) + fcln N\ 
~ S? d (N 2 ,V 2 ,E 2 ) + k\nN\ (22) 

in which the system entropy 5*1 is neglected assuming 
Si <C S 2 . The reduced entropy S 2 ed can be analyzed 
by Taylor's series expansion around N, V, E up to first 
order, that is, 



J 



S? d (N 2l V 2 ,E 2 ) 

S$*(N - Ni,V -V!,E - Ek) 



1 



dS r 2 cd 
dN 2 



S^{N,V,E) 



T 



E,V,N 2 =N 

(—fxNi +E X + V\P) 



E,N,V 2 =V 



~~dE 2 



E 2 =E 



(23) 



with the temperature of the reservoir T. 
The total entropy becomes 



S' = .s'.V-'l.V, V,E) - = (-/iAq +Ex + V X P) + fclnA! 



2 
1 

T 



(-fj,Ni + E 1 + PVi)+C 



(24) 



where C is a constant independent of Ni,Ex,V\. The 
probability <p(Ni,Vx,E\) of the system is proportional 
to the corresponding phase space volume. The phase 
space we are now considering is Tpf 1+ pf 2 , and hence from 
Eq. (JSJ) there are exactly A!/(Ai!A 2 !) copies of the state 
characterized by Nx,Vi,Ei. We have 



piN^VuEi) cx 



AH 



Ai!A 2 ! 



n 1 (N l ,V 1 ,E 1 ) 



1 



cc —^(N^Ei) 

1 fiNi-Eh-PVt 
oc ^ y exp( ) (25) 



under the condition that Ni <C N 2 . 

The probability distribution allows us to formulate the 
Grand Canonical ensemble in which our system inter- 
acts with a particle-energy reservoir. The grand partition 
function 



oo 

E(,,V,T) = EEl 



N r =0 a 



NJ. 



_ e ^N r /kT e -E s /kT 



= J2 e^/ kT Z Nr (V,T) 



(26) 



N r —0 



with the partition function 



(27) 



of the corrected Boltzmann counting due to the factor 
1/AU 
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CONCLUSION 

In summary, although the phase spaces of the system 
of distinguishable particles is different from that of in- 
distinguishable one, their thermodynamics are in fact 
equivalent. It also implies that the corrected Boltzmann 
counting factor can be justified in classical statistical me- 
chanics. 
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